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Abstract 

The spectral principle of Connes and Chamseddine is used as a starting point to define a discrete 
model for Euclidean quantum gravity. Instead of summing over ordinary geometries, we consider 
the sum over generalized geometries where topology, metric and dimension can fluctuate. The 
model describes the geometry of spaces with a countable number n of points, and is related to 
the Gaussian unitary ensemble of Hermitian matrices. We show that this simple model has two 
phases. The expectation value (n), the average number of points in the universe, is finite in one 
phase and diverges in the other. We compute the critical point as well as the critical exponent 
of (n). Moreover, the space-time dimension J is a dynamical observable in our model, and plays 
the role of an order parameter. The computation of (5) is discussed and an upper bound is found, 
(5) < 2. 
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The quest for a quantum theory of gravity is one of the major goals in contemporary re- 
search. Any attempt in this direction involves the understanding of the space-time structure 
at very short distance. It is generally believed that at the Planck scale space-time may not 
be described by a manifold The conventional geometrical setting of general relativity 
seems to be inadequate to describe the non-manifold micro-structures of space-time. Since 
the manifold structure has to appear at some macroscopic limit, it is natural to expect that 
one needs a generalization of ordinary geometry, such as noncommutative geometry (NCG), 

rather than a completely new formalism. 

Q 

The basic idea coming from NCG [2j is that one can describe a Riemannian manifold 
(M, g^) in a purely algebraic way. There is no loss of information if, instead of the data 
(M, g^ u ), one is given a triple (A, TC, D), where A is the C*-algebra C°(M) of smooth func- 
tions on M, TC is the Hilbert space of L 2 -spinors on M, and T> is the Dirac operator acting 
on 7i. From the Gelfand-Naimark theorem it is known that the topological space M can 
be reconstructed from the set A of irreducible representations of C°(M). Metric is also 
encoded, and the geodesic distance can be computed from T>. Here we will consider only 
commutative spectral triples - this is enough to go much beyond ordinary geometry. In 
particular one can treat all Hausdorff topological spaces in this way. Given a pair (M, g^), 
one can promptly construct the corresponding triple (C°(M), L 2 (M), T>). However, not all 
commutative spectral triples, or generalized geometries, come from a pair (M, g^). Never- 
theless one can always associate a Hausdorff space M = A to a commutative spectral triple, 
where A denotes the set of irreducible representations of A. However, the space M may not 
be a manifold. Once we trade the original Riemannian geometry for its corresponding com- 
mutative triple we need a replacement for the Einstein-Hilbert action Seh- The so-called 
spectral action of Chamseddine and Connes [3] is one possible candidate. It depends only 
on the eigenvalues of T> and contains Seh as a dominant term. In this paper however we 
shall use another spectral action. 

The spectral action can be written for any triple, regardless of whether it comes from a 
manifold (M, g^) or not. In the spectral geometry approach it is conceivable to write the 
partition function 

Z = 5>-*M (1) 
were the "sum" is over the set X of all possible commutative spectral triples and S depends 
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on the spectrum of T>. It includes all Hausdorff spaces and therefore all manifolds of all 
dimensions. 

The framework of spectral NCG is appropriate to explore some difficult questions in 
quantum gravity. For example, Rovelli used a simple dynamical model based on a finite 
dimensional spectral triple to construct a quantum theory with a quantized physical distance 
between points 4]. In our model we look at the dimension as a dynamical quantity and 
study its expectation value. 

The exact computation of (1) is a formidable task, not yet accomplished. However, the 
algebraic approach suggests ways of defining discrete approximations to the full theory. For 
instance one may replace the algebra A by a finite dimensional algebra A n . In this approach 
to discretization there is no need to introduce a lattice or simplicial decomposition of the 
underlying space. The approximation of A by a finite algebra works even if the spectral triple 
does not come from a manifold. Thus, it gives us a generalization of ordinary discretizations 

0. 

In this paper we discretize (J^) by sampling the set X with finite commutative spectral 
triples. We will think of it as a useful toy model, which we believe captures some of the 
main features of the full one, Eq. (1). 

The key role played by the eigenvalues of the Dirac (or Laplace) operator in the spectral 
action approach was emphasized in Q|. In our model they are also the natural dynamical 
variables due to the connection with random matrix theory (RMT). An important ingredient 
of the model is that the number of points can fluctuate. Moreover, in our simple model the 
space-time dimension is a dynamical observable and its expectation value can be computable 
from first principles. 

Let us describe the ensemble X C X of geometries we will consider. A point of x G X 
is a commutative spectral triple x = (A, TC, D) where the commutative C*-algebra A has a 
countable spectrum A. We divide X into subspaces X n consisting of triples (A n , TC n , D) such 
that A n has a fixed number n of points. From the Gelfand-Naimark theorem it follows that 
elements of A n are the (possibly infinite) sequences a = (ai, 02, a n ), a,j G C. The Hilbert 
space 7i n is given by vectors v = (v\, ...,v n ) with norm \ \v || 2 = 'Y^ =l vf < 00. The elements 
of A are represented by diagonal matrices a = diag( ) acting on Ti n . Finally, the 

operator D is a n x n self-adjoint matrix. We will sample the space X by Xi,X2, ...,Xn 
and eventually take the limit N 00. 
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Let L be a length scale such that the operator T> given by T> = D/L will be the analogue 
of the Dirac operator. The Chamseddine-Connes action depends on a cutoff function of the 
eigenvalues of D/L. The cutoff function is zero for eigenvalues of D greater than L and one 
otherwise 

a a. 

In other words, the Boltzmann weight in Eq.Q would be one outside a 
compact region in the eigenvalue space, leading to a divergent partition function (see Eq. 
(fTTj) . Let us consider a quadratic action instead: 

^^Trp 2 ), (2) 

where A is the inverse of Planck's length l p , and (3 = (l p /L) 2 . Finally, we define the partition 
function Z N ((3) = J2n=o z n(P) where 

zM = [ [dDje-P^ (3) 



is the partition function restricted to X n , in other words, an integral over all independent 
matrix elements D^, where [dD] is the usual measure for nx n Hermitian matrices |8(. The 
partition function z n (/3) defines the one-matrix Gaussian Unitary Ensemble (GUE) A 
straightforward computation gives z n {(3) = 2^{-^) s ^ . 

The expectation values of an observable O(Dij) restricted to X n and for the entire en- 
semble are 

(0) n ,p = J [dD]Oe-^ D ^ / z n {(3) and (4) 

N 

(O){0) = J2 p (n,(3)(0) n ^, (5) 



n=l 



respectively, where the function P(n, (3) = ^ n( f\p) 1S interpreted as the probability of having 
a universe with n points. The simplest observable in our model is n, the number of points 
in A. By definition, n is constant in X n , therefore (n) n p = n. Thus we get 



E n n2f( 2 , 
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(n)(/3) = ^ • (6) 

The mean (n) ( "average number of points in the universe" ) is not a continuous function of 
(5 at /3 C = 7r/2, signaling the onset of a phase transition. Besides straitforward numerical 
calculation, there are other ways to show that the sum © converges for (3 > (3 C and diverges 
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for (3 < (3 C . Consider, for instance, the approximation of Zn{/3) in the limit N — > oo by the 
first term of a Euler-Maclaurin expansion, 

Z{(3) = lim Y, Z M = / dxe-^ 2+bx + K([3), (7) 
~ Jo 

where ap = ln(/3//? c ), b = In 2/2. It is easily seen that 

(n)(f3) = ^ In Z(0), (8) 

neglecting the remainder 1Z in (|7J) [9J. Eq.((7J) suggests a nice interpretation. For [3 > /3 C 
{(3 < (3 C ) the quadratic term is positive (negative) and the integral converges (diverges). 
The phase transition at (3 = /3 C is triggered by the change of signal of a bilinear term in the 
"fields" . The integral in (J7J) can be solved in the region (3 > (3 C . After expanding ap around 
the critical point, ap = {(3 — (3 C )/ ' (3 C + we compute (n) by means of (JHJ). The result is 
(n){(3) ~ (^-p^j = t1 - Hence, for (3 > (3 C the system is in a "finite" phase, characterized 
by a finite value of (n). As (3 — > (3+ (n) diverges with a (mean field) critical exponent v — 1. 
The r.m.s. deviation of (n) may be computed, with the result An = a/ (n 2 ) — (n) 2 ~ t -1 / 2 . 
However, the relative width of the distribution, An/(n), decreases like r 1 / 2 . 

For f3 < j3 c the relevant universes have (n) = oo and An/ (n) = 0. For a oo- 
dimensional D one can define the dimension 5 of the space A from the eigenvalues of D. 
Let {/io(-D), /ii(-D), •••} be the modules of the eigenvalues (i.e. the singular values) of D 
organized in an increasing order. By the Weyl formula Q], the dimension 5 is related to the 
asymptotic behavior of the eigenvalues for large k: /ifc(-D) ~ ks. By definition 5 = for 
finite dimensional spectral triples. We can argue that (5) is of the form 

(0 if/3>/3 c . 

This follows from the fact that for (3 > (3 C the probability P(n,/3) is localized around some 
finite n. Hence (5) works as an order parameter. The value (3 C = n/2 separates (5) = from 
the rest. 

In order to study the dimension we need to consider the spectral ^-function 

n 

C(z)= hm = (pi"*), (10) 

k=0 
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where D is an cx)-dimensional matrix (/io > 0). The relation between the dimension and ((z) 
has been discussed in For large enough values of a = Re (2), Tr(|.D|~ 2 ) is well defined. 
One says that D has dimension spectrum Sd if a discrete subset Sd = {s\, S2, ■■■} C C exists, 
such that ((z) can be holomorphically extended to C/Sd. This definition is consistent with 
the Weyl formula. The set Sd has more than a single point when for example the geometry 
is the union of pieces of different dimensions jll|. In what follows we will look at an upper 
bound for the dimension: It may happens that Tr(|-D|~") = for large enough a, whereas 
for small values of a, Tr(\D\~ a ) = 00. Eventually, there is a value of a (say, a c ) for which 
Tr(|/}| -Qc ) is finite and non-zero. The upper bound for the dimension will be 5 = a c . 

In order to estimate {5) by means of (|1(J|). we rewrite (jSJ) and (@J) as integrals over the 
eigenvalues A^ of D. The procedure is well-known j8|, and leads to (C n = n ( 2 ) / YYk=i 

/oo 
[d n \]A 2 (\ k )e-^^ x * = C n * n ,p , (11) 
-00 



{0{\)) n ,p = / [d n X]0(Xi) \ n M A 2 (A fc )e-^ A * = / [^©(A^P^^Afc) , 

J-00 [ 7T 2 llfc=l fc! J ./-oo 

(12) 

where A(A^) = Yli<j(^j~ ^i) ^ s ^ e Vandermonde determinant (Jastrow factor), and [d n \] = 

In RMT, ^n.pil) is interpreted as the positional partition function of an ensemble of 
equal charged particles (with positions given by Aj) in two dimensions, moving along an 
infinite line, in thermodynamic equilibrium at temperature 7 - the so-called "Dyson gas" 
[lfl| . Then, V n ,p{X\, A n ) defined in (fT2*|) is the probability of finding one particle at Ai, 
one at A 2 , etc. The value of ^^(7) is known from the Selberg's integral. 

In the region (3 < f3 c the partition function Z{f3) is dominated by 00— dimensional ma- 
trices. Thus, one may try to select the 00— dimensional matrices out of the whole ensemble, 
and then compute the mean of the ^-functions following (JTUJl . However, from the standpoint 
of our statistical approach this procedure does not seems natural since the sum over n is 
a key ingredient in the whole construction. Hence, we look for a quantity related to the 
^-function that captures the statistical nature of our model. Let us compute the mean value 

(Tr K \D\- a ) n>f3 = (j2 \X k \- a e{\\ k \-K)\ . (13) 

\fc = l / n R 



Bearing in mind (|10jh we want to study (Tr K \D\~ a ) n ^ only in the region of eigenvalues where 
it is a decreasing function of a. Hence, we introduced an "infrared" cutoff k, since there is 
a non-zero probability to find a configuration in the volume ^ \Xk\ < k around the origin. 
Using the symmetry of the integrand in (|13|) under permutations of the position indexes, we 
arrive at 

(Tr K |L>|~ a ) ni/3 = f dX n \X n \- a a n (y^X n ), (14) 

J OUt,K 

where we have used the definition of the spectral density a n (A) = (Y^ S(X — Xk)) n - Besides, 
J o = J_^ + ■ The computation of (|T3|l was reduced to a one-particle problem. It is 
known that cr n (A) = J2 k <Pk(^)> 4>k{ x ) being the Weber-Hermite functions ^]. In the large 
n-limit, cr n (X) converges to a nonrandom function, Wigner's "semi-circle law", 



Mv^ah^^i-I^, (is) 

for |A| < 1Zn t p = \f2nf (3, and zero otherwise. The average eigenvalue density for Gaussian 
ensembles of Hermitian matrices with different values of n differ only by a change of scale, 
for large enough n (since TZ n ^ ~ V^)- We propose the following quantity in order to extract 
information on the (^-function ()10|): 

N 

(Tr Kn jDD([3) = lim £ P(n, /?)(Tr K JI}|- Q ) n>/3 . (16) 

n=l 

Hence, we are sampling the partial traces of the ensemble of oo— dimensional Hermitian 
matrices by the total trace of finite dimensional matrices in the complete Gaussian ensemble 
of Hermitian matrices. Now, we study the asymptotic properties of (Tr Kn/3 \D\~ a ) n ^ as a 
function of a and n for n ~ N large. In this case, one may use the semi-circle law in (fT4l) . 

(Tr Kr jD\~ a ) rh(3 « -V2n/3 / rfA n |A n |- Q a n (v^A n ). (17) 

In (|T7j) we are neglecting the contribution from the "tail" of <r n (y/j3X), outside of the semi- 
circle radius TZ n ,/3- There is an exponential decrease in the tail, so that the total number 
of particles (or eigenvalues) in this region is of order 2 [12]. This finite-size correction is 
not relevant to our asymptotic analysis. Besides, the cutoff k is a function of n and (3: 
K n fi = e^2n/ (3. The mean particle (or level) spacing is s(A, (3) ~ a^ 1 (X,(3), so that the 
mean spacing is almost uniform near the origin. Near the endpoints (edges) the spacing is 
highly non-uniform. This is the asymptotic region we are interested in. Our procedure is to 



select a slice of size A n) p )K = H n> p — n n ,p = (1 — ^)T^n,/3 near the edge, so that the relative 
size of the slice, A. n ^ K /TZ nt /3, does not depend on n and (3. This choice ensures that we are 
treating matrices of different sizes n on the same footing. From (|17j) and (|14|) we obtain, 
after some manipulations, 

<Tr KjW3 \D\-"U « I (2n) 1 -f 0f jf * ^ v 7 ^- (18) 

The asymptotic behavior of (Tr Kn(3 |.D| _Q! ) nii g does not depend in an essential way on the 
particular choice of e, as long as we keep e ^ 0. 

Now we use the asymptotic formula (j!8)l in (jlfij) and search for the value a c for which, 
as iV — > oo and /3 — > /5 C , (Tr Kn(3 |-D| _a ) diverges (converges to zero) if a < a c (a > a c ), 
with (Tr Kn/3 \D\~ ac ) finite and non-zero. This gives an upper bound for the dimension of the 
"condensed" manifold in the infinite phase ( (3 < /3 C ), which is {5) < a c . We obtain: 

N 

(Tr Knt0 \D\- a )((3) ~ lim ^ P{n^)n 1 ^ . (19) 

n=l 

In the finite phase (/3 > (3 C ) the sum in (flDj) converges for a > 0. We conclude that a c = 
(i.e. (5) = 0) for /3 > f3 c , as expected. From the behavior of P(n, (3) in the infinite phase it 
follows that the convergence of the sum in ()19j) is dictated by the behavior of T n ^ a = n l ~^ 
in the limit n ~ N —>■ oo. For (3 < /3 C we get r n . a — > oo if a < 2, and T n ^ a — * if a > 2. 
For a = 2 it turns out that (Tr Kn/3 \D\~ 2 )(j3) ~ 1. Therefore, we obtain the upper bound 
(5) < 2. 

Notice that we do not have a definition of the dimension as an operator (as suggested by 
Q), from which it will be possible to compute its average in the ensemble of all Hermitian 
matrices, like we have do for n in Instead, our guide is the operational definition 

encapsulated in In order to go beyond the upper bound computed above we need to 

find a suitable observable which reduces to the classical dimension (as given by the Weyl 
formula) in the limit n — > oo. This eventually may lead to a numerical computation of (5). 
Besides, one sees that the reason for the upper bound (5) < 2 lies in the semi-circle law 
and its leading \/A behavior near the edge. It is known from 2D models of discretized pure 
quantum gravity that some special matrix polynomial potentials V(D) = akTrD 2k 

may lead, in a suitable scaling limit, to a behavior near the edges different from the square 
root. Thus, a possible way to obtain a bound of higher dimension would be to include 
higher polynomial interactions, or modify the quadratic one in (J2J) including more (internal) 
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symmetries besides the unitary one. These and other questions are under study, and will be 
reported elsewhere. 

To conclude, in this Letter we proposed a discrete model for quantum gravity based 
on the spectral principle. The model is connected with the GUE of Hermitian matrices, 
and contains the mean number of points, (n), and the dimension of the space-time, (5), as 
dynamical observables. We have shown that the model has two phases: a finite phase with 
a finite value of (n) and (5) = 0, and an infinite phase with a diverging (n) and a finite 
(5) 7^ 0. The critical point was computed, (3 C = tt/2, as well as the critical exponent of (n). 
Moreover, an upper bound for the order parameter (5) was found, (5) < 2. 
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